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SUMMARY 


nr** 


Equations which are easily adaptable to machine computation are presented 
for the preliminary prediction of the supersonic and hypersonic stability 
characteristics of bodies of revolution with an arbitrary center -of -gravity 
location. The equations have been programed and an example of their application 
to a blunt reentry configuration is presented. 


INTRODUCTION 


The Newtonian theory has proven to be a useful tool for the preliminary 
prediction of supersonic and hypersonic stability characteristics. A general 
background in the application of the theory may be obtained from references 1 
to 4. The machine computation of the Newtonian static stability characteris- 
tics of bluff entry bodies with various heat-shield curvatures, afterbody 
angles, and corner edge radii is in agreement with experiments on Apollo type 
reentry configurations. (For example, see ref. 5.) Since entry lift-to-drag 
ratio may be obtained by offsetting the center of gravity, the program was 
extended to the computation of both static and dynamic coefficients for bodies 
of revolution with an offset center of gravity at angles of attack, angles of 
sideslip, and angles of roll. The purpose of this paper is to present equations 
which are readily machine adaptable for the computation of the Newtonian static 
and dynamic aerodynamic coefficients for a body of revolution with an offset 
center of gravity at angles of attack, sideslip, or roll, or any combination 
of these angles. Q 
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defined algebraic expressions 
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axial -force coefficient. 
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yawing-moment coefficient. 
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roll, yaw, and pitch damping coefficients 
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normal force 

unit vectors for tody coordinates 
lift -to -drag ratio 
Mach number 
rolling moment 

pitching moment 

yawing moment 

unit vector normal to surface, positive inward 

surface vectors in roll, pitch, and yaw planes 
roll, pitch, and yaw rates about body-axis system 

free-stream dynamic pressure 
body radius 

relative velocity, normal to surface 

free-stream velocity 

body coordinates 
distance along body sixes 
center-of -gravity coordinates 

x 

lower and maximum integration limits for — 

angle of attack 
angle of sideslip 
surface slope from "body axis , 6 

cylindrical coordinate angle 
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0 j 

flow-see boundary limits 

®a 

body afterbody angle 

P 

radius vector for cylindrical coordinates 

0 

roll angle 

Subscripts 

(see fig. k): 

a, v 

afterbody - vertex juncture 

c 

corner edge 

c,a 

corner edge - afterbody juncture 

n 

nose 

n,c 

nose - corner edge juncture 

V 

vertex 

vr 

vertex radius 


DEVELOPMENT OF EQUATIONS 

The Conventions and Basic Relations 

The vector directions used in this paper are found in the diagram of 
figure 1. The respective yaw, pitch, and roll rate vectors (see fig. 2) are 


r = - rk (la) 

(lb) 

P = -Pi* (1c) 

The surface vectors normal to their respective rotation axes (fig. 2 ) are 

R" = (x-x 0 )i* - p cos Q 3* ( 2a ) 

( X “ X 0 ) I> + (p sin e - z 0 )? (2b) 

p* *= - P cos e3* + (p sin e - z Q )k* ( 2c ) 
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The inwardly directed unit vector , normal to the surface, is 


— = sin 5? + cos 5 cos 0j* - cos 5 sin 0k 
|n| 

The unit free-stream velocity vector from figure 5 is 


( 5 ) 


v ^ 

= cos a cos 3? - (sin a cos 3 sin 0 + sin 3 cos 0)j 

I oo| 

+ (sin a cos 3 cos $ - sin 3 sin 0)k 


00 


The relative velocity, normal to the surface, may he expressed as 

-» — > 

n -» 7^ n 

|n| |n| q X Q * |n| 


V, T = t-,-P xP *y^ x ^N- rXR *N 
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is 


The local surface pressure coefficient, according to the Newtonian theory, 


C P ■ 2 (v 
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= 2 (A cos 0 - B sin 0 + C)^ 
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The coefficients may he expressed as 
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Flow-See Boundary 


The flow-see boundary is given by the integration limits 

established at each -7 station. 
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0 ^ and 9,_> 
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The lower limit of the integration with respect to ^ is given by the condition 
where 5 - - cos (cos a cos (3 )J =0. 

For values of ^ not included in the above two paragraphs ^ 0^ and 0^ 

are found by setting the pressure coefficient equal to zero. 



7 


Coefficients 


Substitution of the pressure coefficient (eq. 6) into the coefficient 
equations (eqs. (7a) to (7*)) enables an integration of the coefficients with 
respect to the coordinate 0. The resulting coefficients expressed below 
require only one integration with respect to the axial coordinate of the body 
selected. For other then simple bodies the coefficients may be obtained by a 
machine integration for each combination of a, p, p, q, and r. 
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TYPICAL APPLICATIONS OF EQUATIONS 


The body force, moment, and damping coefficient equations have ^een 
programed on an IBM 7 O 9 O computer. The typical blunt reentry configuration 
geometry, as presented in figure 4, has been programed into the body coeffi- 
cient expressions. The body equations for the configuration shown in figure 


are: 


For 
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n, c 


d 


then 


and 



(11a) 


(lib) 


11 
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A computed case for a typical blunt reentry configuration with — = 1 . 2 , 

R R 

— = 0. 05, 0 = 35°, and -- = 0.1 is compared with theory in figure 5. 

d ’ a a 

The constants for the configuration shown in figure 4 for this case are: 
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DISCUSSION 


Typical Application 

Results . - Generally good agreement with experiment was obtained for both 
the static and dynamic coefficients , as shown in figure 5 . The experimental 
results for the static coefficients were extracted from unpublished data 
obtained in the 21-inch hypersonic tunnel of Jet Propulsion Laboratory, 
California Institute of Technology. The experimental results for the dynamic 
coefficients were taken from unpublished data obtained in the Langley Unitary 
Plan wind tunnel. 

Approach to indeterminateness . - Some of the expressions, as presented, 
are invalid where the body surface is normal to the body axis of revolution 

p 

because of the appearance of the indeterminate form — tan 5 . These points may 
be singularly evaluated or values may be obtained by extrapolation. For the 
example configuration as — -» Oh: -» 0 , tan 5 the value of ^ tan 5 is 
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B 

which may he substituted into the expressions. However, for the compu- 
tations in figure 5, when applicable the integration was started and stopped 

X 

0. 0001 — interior of the points where the hod y slope is normal to the hod y of 
revolution; hence, avoiding the indeterminate points. 


General Application 

Choice of stagnation point coefficient . - All the coefficients computed 
from the presented equations may he modified to correspond to the actual 
stagnation point pressure coefficient hy multiplying the computed coefficients 
hy the ratio of the actual stagnation point pressure coefficient to the 
Newtonian value (2.0). 

Extension of application . - It may he noted that the expressions in this 
paper are only as valid as the assumptions of Newtonian or modified Newtonian 
Flow Theory. The accuracy of this theory depends on the particular con- 
figuration under consideration. 


CONCLUDING REMARKS 


Equations which are easily adaptable to machine computation are presented 
for the preliminary prediction of the supersonic and hypersonic stability 
characteristics of bodies of revolution. The equations have been programed 
and an example of their application to a blunt entry configuration is given. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January JO, 1963 
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Figure 4 . - Configuration setup for programing into tody coefficient 

expressions. 
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